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Abstract 

In this paper we present a hybrid model of k-essence and chameleon, named as k-chameleon. In 
this model, due to the chameleon mechanism, the directly strong coupling between the k-chameleon 
field and matters (cold dark matters and baryons) is allowed. In the radiation dominated epoch, the 
interaction between the k-chameleon field and background matters can be neglected, the behavior 
of the k-chameleon therefore is the same as that of the ordinary k-essence. After the onset of matter 
domination, the strong coupling between the k-chameleon and matters dramatically changes the 
result of the ordinary k-essence. We find that during the matter-dominated epoch, only two kinds 
of attractors may exist: one is the familiar K attractor and the other is a completely new, dubbed 
C attractor. Once the universe is attracted into the C attractor, the fraction energy densities of 
the k-chameleon fi^ and dust matter $1™ are fixed and comparable, and the universe will undergo a 
power-law accelerated expansion. One can adjust the model so that the K attractor do not appear. 
Thus, the k-chameleon model provides a natural solution to the cosmological coincidence problem. 

PACS numbers: 98.80.-k, 98.80.Cq, 04.80.Cc 



* Email address: haowei@itp.ac.cn 
t Email address: cairgOitp. ac.cn 



I. INTRODUCTION 



There are now a lot of cosmological observations, such as SNe la HHI3, WMAP SDSS etc.. AU 
suggest that the universe is spatiaUy flat, and consists of approximately 70% dark energy with negative 
pressure, 30% dust matter (cold dark matter plus baryon), and negligible radiation, and that the universe 
is undergoing an accelerated expansion. To understand the nature of the dark energy remains as one of 
biggest challenges to theorists and cosmologists . The simplest candidate of the dark energy is a tiny 
positive cosmological constant. However, it is difficult to understand why the cosmological constant is 
about 120 orders of magnitude smaller than its natural expectation, namely the Planck energy density. 
This is the so-called cosmological constant problem. Another puzzle of the dark energy is the cosmological 
coincidence problem, i.e. why are the dark energy density and the dust matter energy density comparable 
now? and why does the universe begin the accelerated expansion just only near recent! 

In order to have an interpretation to the accelerated expansion of the universe, many alternatives to 
the cosmological constant have been proposed. One of interesting scenarios is the so-called quintessence 
model 0. The quintessence is a slowly varying scalar field with a canonical kinetic energy term. With 
the evolution of the universe, the scalar field slowly rolls down its potential. A class of tracker solutions of 
quintessence 0,0 is found in order to solve the cosmological coincidence problem. As is shown, however, 
the quintessence model still needs some fine-tuning in order for the quintessence component to overtake 
the matter density at the present epoch (for example, see Refs. 0112). Motivated by the k-infiation [l^, 
in which a scalar field with non-canonical kinetic energy terms acts as the inflaton, the so-called k- 
essence j^, 0, is introduced to the coincidence problem. In this well-known model, by the help 

of non- linear kinetic energy terms, a dynamical solution to the cosmological coincidence problem without 
fine-tuning is possible [ill [l2|. In fact, k-essence is based on the idea of a dynamical attractor solution 
which makes it act as a cosmological constant only at the onset of matter domination. Consequently, 
k-essence overtakes the matter energy density and makes the universe start with accelerated expansion 
just near recently. It is worth noting that to achieve a later-time acceleration attractor, one needs to 
design the Lagrangian of the model so that r'^{yd) > 1 in order to avoid the dust attractor jl3|. After 
all, the quintessence and k-essence fields are very light scalar fields. Such light fields may mediate a 
long-range force, and therefore are subject to tight constraints from the searches of the fifth force 
and the tests of the equivalence principle (EP) [l^ . 

On the other hand, the coupling between scalar field and matters has been studied for some years (for 
example, see Refs. [SlllllliEillST', 281). Recently, a novel scenario named chameleon l2(J, 121], 
|2^ ,23s, _24] ha-s been proposed (see also 29]). In this scenario, the scalar field can be directly coupled 
to matters (cold dark matters or baryons) with gravitational strength, in harmony with general expects 
from string theory, while this strong coupling can escape from the local tests of EP violations and fifth 
force searches. The basic idea of the chameleon scenario is that the scalar field acquires a mass which 
depends on the ambient matter density (so the name chameleon). While nearly massless in the cosmos 
where the matter density is tiny, the chameleon mass is of order of an inverse millimeter on the earth 
where the matter density is high, which is sufficient to evade the tightest constraints from the tests of 
EP violations and fifth force searches. 

It is interesting to wonder what will happen when the k-essence is strongly coupled to matters through 
the chameleon mechanism. May the virtues of k-essence and chameleon join together and the shortcomings 
be avoided? The answer is yes. In this paper, we will combine the k-essence with chameleon and present 
a so-called k-chameleon model, which of course is a hybrid of the k-essence and chameleon. Through 
the chameleon mechanism, the directly strong coupling between k-chameleon and matters (cold dark 
matters and baryons) is allowed. We study the cosmological evolution of k-chameleon and find that the 
k-chameleon model can provide a natural solution to the cosmological coincidence problem. 

In the k-chameleon model, during the radiation dominated epoch, the interaction between the k- 
chameleon field and ambient matters can be negligible. Therefore the behavior of the k-chameleon is 
the completely same as that of the ordinary k-essence without the interaction between the scalar field 
and background matters. As a result, three kinds of attractors, namely R, K, S (following the notations 
of the k-essence model 0,0) may exist. The radiation tracker, i.e. the R attractor, has the largest 
basin of attraction on the whole phase plane so that most initial conditions join onto it and then makes 
this scenario became insensitive to initial conditions. However, after the onset of matter domination, 
the strong coupling between k-chameleon and matters dramatically changes the result for the ordinary 
k-essence. In the matter-dominated epoch, the D and S attractors (which may exist in the ordinary 
k-essence model) are physically forbidden due to the strong coupling between k-chameleon and matters. 
Note that unlike the ordinary k-essence model, the disappearance of D and S attractors naturally occurs 
in the k-chamclcon model, need not any artificial design of the Lagrangian. Actually, during the matter 
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dominated epoch, only two kinds of attractors may exist: one is the famiUar K attractor and the other 
is a completely new one named as C attractor. The new attractor C has some desirable features which 
may provide a promising solution to the cosmological coincidence problem. 

Once the universe is attracted into the C attractor, the fraction energy densities of the k-chameleon 
and the matters flm are fixed and they are comparable. Further, many parameters, such as the 
parameter of equation-of-state of the k-chameleon field and its kinetic energy term X , are also fixed. 
And the universe will undergo a power-law accelerated expansion forever. In this sense the k-chameleon 
model gives a natural solution to the cosmological coincidence problem. On the other hand, note that if 
the kinetic energy term X of the k-chameleon is fixed at a somewhat small value (equivalently y = 1/ Vx 
is large), the k-chameleon can be treated as a canonical chameleon approximately. Therefore, we cannot 
detect it from the tests of EP violation and fifth force searches on the earth and in the solar system today, 
although it is strongly coupled to background matters. 

This paper is organized as follows: In Sec. ^ a brief review of the chameleon mechanism is given. 
In Sec. mil we present our k-chameleon model and illustrate how the directly strong coupling between 
k-chameleon and matters (cold dark matters and baryons) is allowed while it cannot be detected from 
the tests of EP violation and fifth force searches on the earth and in the solar system. In Sec. IIVI the 
cosmological evolution of k-chameleon is studied and the result shows that the k-chameleon model may 
provide a promising solution to the cosmological coincidence problem. A brief conclusion will be given 
in Sec. E| 

We use the units h — c — 1 throughout this paper. Mpi = {SttG)^^^^ is the reduced Planck mass. We 
adopt the metric convention as (H ). 



II. A BRIEF REVIEW OF THE CHAMELEON MECHANISM 



Following Refs. |2rll2lLl22| . consider a canonical chameleon scalar field (j) governed by the action 



S = / d^x^ 



1 



+ / d4xA„(V'«,5^)), (1) 



where g is the determinant of the metric 5^^, TZ is the Ricci scalar and ipm are various matter fields 
labeled by i. In the chameleon mechanism, the scalar field is supposed to directly interact with matters 

through a conformal coupling. In other words, each matter field ipm couples to a metric g^l, which is 
related to the Einstein-frame metric 9^,^ by the rescaling 



9U = e'^'^/''^'9,., (2) 

(i) 

where Pi are dimensionless constants. Moreover, the different 'ipm fields are assumed not to interact with 
each other for simplicity. From the action Eq. the equation of motion for </> is 

=-v,-y: -^^''"''^'9^1^, (3) 



where 



V^(P = g^^V^V,^ = -^d^ [^g^^d,^] 



and T^l; — {2/ \/^g^)SCm/Sg'^^^ is the stress-energy tensor density for the i-th form of matter. V^^ 

denotes the derivative of V with respect to (p. For non-relativistic dust-like matter, g^^^Tji^J = pi, where 
Pi is the energy density. Defined in this way, however, pi is not conserved in Einstein frame. Instead, it is 
more convenient to define a matter density p„ii = Pi e^'^''^/^^p' which is independent of cp and is conserved 
in Einstein frame. Thus, Eq. ^ can be recast as 

= -V,-Y: ^Pm. e^^^/^- = -Vf. (4) 
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Note that the dynamics of (f) is not governed solely by V{4)), but rather by an effective potential 



i 

which depends explicitly on the matter density pmi ■ The key ingredient to achieve a successful chameleon 
model is that the effective potential Veff{4') has a minimum even when V{(l)) is monotonic. In fact, if V{4i) 
is monotonically decreasing and (3i > or, equivalently, V{<j}) is monotonically increasing and j3i < 0, the 
effective potential T4/ / (0) has a minimum 0mm satisfying 

= V^^rmn) + V Prm e^^'^^-/^'^' = 0. (6) 

^—^ Mpi 

Meanwhile, the mass of small fluctuations about the minimum (pmin is 

miff ^ V;//{cl,^.n) = F^40™„) + -^prm e^^^--/''^^. (7) 



In other words, the originally massless scalar field acquires a mass which depends on the local matter 
density. The denser the environment, the more massive the chameleon is. Actually, while the coupling 
constants f3i can be of order unity as the natural expectations from string theory, it is still possible for 
the mass of chameleon, i.e. We//, to be sufficiently large on the earth to evade current constraints on 
EP violation and fifth force. On the other hand, through the so-called "thin-shell" effect, the chameleon- 
mediated force between two large objects, such as the earth and the sun, is much suppressed, which 
thereby ensures that solar system tests of gravity are satisfied. For more details, see the ori gina l papers [2(1 
121^ .22>. ,23] . The quantum stability analysis of the chameleon model is presented in Ref. ]3| • 

It is worth noting that, in most existing canonical chameleon models, the potentials V{(j)) are assumed 
to be of the runaway form, namely it is monotonically decreasing and satisfies 

hm V = 0, lim ^ = 0, Um = ... 



as well as 



lim V = oo, lim — ^ = cxd, lim 



+0 V 
Actually, the fiducial potentials are chosen to be 

V(0) = i—j and = M"* exp(M"/0") (8) 

in Ref. [2(1 12^ and Ref. |23|, respectively. However, for the potentials given in Eq. ©, to achieve a 
successful chameleon model, the mass scale M has to statisfy 

M < 10-^ eV, (9) 

which is about 30 orders of magnitude smaller than its natural expectation, namely the Planck mass. 
Fortunately, Ref. [2^ shows that a chameleon model with a non-runaway form potential 

can be a successful example even when the parameter <^ is of order unity. In Sec. IIII CI of the present 
paper, a completely new canonical chameleon with non-runaway potential and without any fine-tuning 
like in Eq. will be presented. This new canonical chameleon model is another successful example. 



III. K-CHAMELEON MODEL 



So far, all chameleon models existing in the literature is of the form of quintessence-like, namely the 
kinetic energy term of the scalar field is canonical one. As is well-known, non-linear kinetic energy terms 
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naturally appear in many models unifying gravity with other particle forces, including supergravity and 
superstring theory. For many years, the contributions of these hi ghe r order terms have been ignored 
for the reasons of simplicity. The example of k-essence [ll[ [lj| demonstrates that the effects of 

non-linear dynamics can be dramatic. Here, motivated by k-essence, we put the chameleon and k-essence 
together and present a k-chameleon model. Namely, we consider a scalar field with non-linear kinetic 
terms and the scalar field is strongly coupled to matters. 



Our starting point is the action 



(Tx 



Af2 



A. Setup 



(11) 



where g is the determinant of the metric 5^,^, TZ is the Ricci scalar and ipm are various matter fields 
labeled by i. The scalar field (j) interacts directly with matters through a conformal coupling. Explicitly, 
each mat1 
rescaling 



each matter field tpm couples to a metric glil which is related to the Einstein-frame metric g^i, by the 



g<^^AIi^)g 



(12) 



where Ai{<j)) is a function of the k-chameleon field <j>. Note that Eqs. Hll|l and H12|l are of the general 

form arising from string theory, supergravity and Brans-Dicke theory. Moreover, the different fields tpm 
are assumed not to interact with each other for simplicity. The kinetic energy term is defined by 



X = -g^^d^c^d,^ 



(13) 



Note that if one takes p((/>, X) = X — V{4)) and Ai{(f)) = exp(/3i0/Mp;), this case then reduces to the 
canonical chameleon model presented in Sec. UTI 

Varying the action Eq. (|ll|l with respect to (j) yields the equation of motion for the k-chamclcon field 



-9 

where p.^ denotes the derivative of p with respect to X , and 

2 S£m. 



is the stress-energy tensor density for the i-th form of matter, and 

"*^^) = ^^- 

For the i-th non-rclativistic dust-like matter, the energy density in Einstein frame 

2 SCr. 



p,=T^^ = g^^- 



f^5g^^'■ 



where 



2 SCm 



(14) 



(15) 



(16) 



(17) 



(18) 



is the energy density in the matter frame. However, pi is not conserved in Einstein frame. Instead, it is 
more convenient to define matter density 



(19) 
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which is independent of (j) and conserved in Einstein frame. Then Eq. (|14|l can be rewritten as 

-^d,, [V^p,^ g^-d^cb] = P,. - E "»('^) P^- (20) 
In addition, form the action Eq. we have the pressure and energy density of the k-chameleon field 

p^^p{(j),X), p^ = 2Xp^^ - p , (21) 
respectively. Consider a flat FRW universe, whose metric is 

ds^ ^ dt^ - a^{t)dyi^, (22) 
where a is the scale factor. If the scalar field (/) is spatially homogeneous, one then has 

X = (23) 

where a dot denotes the derivative with respect to the cosmic time t. Furthermore, by using Eq. H21|l . 
the equation of motion for the k-chameleon field 0, namely Eq. (|20|) . can be recast as 

p^ + iH{p^^p^) = -^a^{(j)) pi(j), (24) 

i 

where H = a/a is the Hubble parameter. From the total energy conservation equation 

ptot + iH{ptot + Plot) = 0, (25) 
where the total pressure and energy density are 

Ptot =p^+Pr, and Ptot = P4, + '^Pi + Pr, (26) 

i 

respectively, we have 

Pi + 3Hpi=ai{(l))pi<j), (27) 

and 

Pr + AHpr = 0, (28) 

where Pr = Pr/i and pr are the pressure and energy density of radiation respectively. Note that there 
is no coupling between the scalar field and radiation because the trace of the stress-energy tensor of 
radiation vanishes. Finally we write down the Friedmann equation 

2,H'^Mli = Ptot =P^+ Pm + Pr (29) 

where Pm = X] Pi sum of the energy density of all matter components. 

i 

B. Master equations of k-chameleon 

In this paper, following k-essence jlllll2j |. we only consider a factorizable Lagrangian of the form 

p{^,X)^K{^)D{X), (30) 

where we assume K{(l)) > 0. From Eq. H21() . we have 

p^^p{(^,X)^Ki4>)D{X), 
p^ = 2Xp,, -p^ K{cf>) [2XD,AX) - D{X)] = K{4,)E{X). (31) 
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The parameter of equation-of-state is 



P4> 



D 

E 



D 



2XD ^ - D 



According to the definition in jlflj, the sound speed is 
Substituting Eqs. lp?T|l and (|^ into Eq. it^ . we obtain 



dx 

dN 



E 

eZ 



3(1 + w^) + G 



K H 



V2X 
HKE 



(32) 



(33) 



(34) 



where = In a is the so-cahed e-folding time, a is the sign of ((>. 

It is convenient to re-express D as D = g{y)/y and to view it as a function of the new variable 
y = l/\/X. Eqs. (EOll-ljSSI) then become 



W4, 



Kg 

1 

y 
9 

g'y' 



Prp = -Kg', 



g,fy2 



(35) 

(36) 



where a prime denotes the derivative with respect to y. Taking into account Eqs. (|29|l . H35|) and H36|) . 
one can recast Eq. ()34|l in terms of the new variable as 



dy_ _5 w^{y) - 1 
dN ~ 2 s'{y) 



V Ptot 



where 



s{y) 



5g' 



1/2 



yi^ + w^) 



g y 



-9' 



(37) 



(38) 



Note that the requirements of positivity of the energy density, > 0, and stability of the k-chameleon 
background, > 0, imply 



-g' > 0, 



g" > 0. 



(39) 



These conditions indicate that g should be a decreasing convex function of y = 1/VX. A sample of the 
function g{y) is plotted in Fig. 1. 

Since we are attempting to study the cosmological evolution of k-chameleon and try to find out its 
attractors, we impose the condition that the coefficients of the last two terms in Eq. 137|l to be constants 
for simplicity. Thus, if(0), ai{(j)) and A.i{<j)) should be the form 



1^' 



a,,((/)) 



f3^ 



A.(0) 



M J 



(40) 



where M is a constant mass scale and Pi are dimensionless positive constants. Furthermore, although Pi 
may be different for different matter species, we take a same value P for all Pi for simplicity [«(</>), A((f>) 
for all ai{(j)), Ai{(j)) accordingly]. In fact, it is straightforward to generalize it to the generic case with 
different Pi . In addition, without loss of generality, we restrict ourselves to the most interesting case of 
positive (j), namely a = +1. Under these simplifications, Eq. H37(l becomes 



dy_ 3 w^{y) - 1 
dN ^ 2 r'{y) 



riy) 



pn^ 
2 ./Th, 



(41) 



where r(y) = Ms{y) is a dimensionless function of y. It is worth noting that if M = ^/SMpi, r{y) then 
is completely the same as that of the k-essence case 0, . The fraction energy densities fi^ = p^/ ptot 
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and — Pm/ptot, where p„i — pi is the sum of the energy density of all matter components. At the 

i 

same time, from Eqs. ()24f) - (|28|l . we reach 

flm4> + S^ct>iwtot - W^), (42) 



dN H 



dN H 



ilrn'i' + i^mWtot, (43) 



where 



Wtot = = w^n^ + l^^r, (44) 
Ptot o 

and Qr = Pr/ Ptot is the fraction energy density of radiation. 

Thus we have obtained the master equations governing the whole system. Compared to the ordinary k- 
essence model |lj||12j, we find that some new terms which describe the coupling between k-chameleon and 
matters enter these equations. In addition, let us mention that in our k-chameleon model, the c oupling 
function A{(j)) is of a form of power-law [see Eq. (|4()|l ] . while in the ordinary chameleon models [2(ll2lll2^ 
l2^ 0| or other coupled models of scalar field to matters in the literature jl^ , the coupling function 
is of either the exponential type [H HO, IH HI [H HI] or linear type pLOj . 



C. Coupling constant /3 and mass scale M 

In this subsection we will discuss the constraints on the coupling constant f3 and the mass scale M 
imposed by the fifth force experiment and the tests of the equivalence principle (EP) in order 
to obtain a successful chameleon mechanism. 

Note that all the searches of the fifth force and EP violation have been performed only at the present 
epoch on the earth or in the solar system. This point is very important. As we will see, in the k- 
chameleon model, the kinetic energy term X of the k-chameleon field could be fixed at a somewhat 
small value (equivalently y = is large) at the current accelerated expansion epoch. In this case, 

the non-linear kinetic energy terms can be neglected and the k-chameleon can be treated as a canonical 
chameleon approximately. 

In this case, D{X) in the Lagrangian p(0, X) = K{(j))D{X) can be approximated to 

D{X) C:i CiX ^ C2, (45) 

where ci is a dimensionless positive constant and C2 is a positive constant with dimension of energy 
density. To change the Lagrangian to a canonical form, we make a redefinition of the field variable. 
Introducing a new scalar field 

(j) 

4>new = -\/c\M In (46) 



satisfying 



^neu) = -j^g^" d^j,(j)newdv4>new = CiK{(j))X, 



where K{(f)) = has been considered. Accordingly, the conformal coupling A(0) = [(p/M)^ is 

transformed to 



A{(t)new) = exp 



P4>ne 



(47) 



and the Lagrangian becomes the canonical one, 
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where the potential 



V(4>new) = C2exp 

The equation of motion Eq. 120|) then becomes 



2(j)„ 



where the effective potential 



and the matter energy density 



Veff{4>new) = V {(pnew) + A{4>new) P\ 



new 
'rat ' 



^new ^3 / 



Pmt 



(49) 

(50) 
(51) 
(52) 



which is independent of ipnew and is conserved in Einstein frame, and pi is defined in Eq. IjlSI) . Obviously, 



because V 



is always larger than zero, the effective potential has a minimum at 



rmin 
^new 



2C2 



2 + (3 



satisfying V^J^ (</'™e™) = 0. Thus, the mass of small fluctuations about the minimum 0™™ is 



V 



eff 



ini'in 
^new 



)• 



After some algebra, we get 



Cl 



1/2 



M 




^new 



-l/(2+/3) 



,/?2 + 2/3y \M^i) \ 2c2 
It is easy to see that, if one takes the natural expectations of the constants as 

Cl ~ ^ - 0(1), C2 ~ C(M^,) and M - 0{M^ 



m: 



(53) 



(54) 



(55) 



(56) 



"^new is indeed a small quantity for any reasonable matter density. For instance, the atmosphere has 



mean density Patm ~ 10 g/cm'^. Substituting into Eq. (|55|l and assuming Eq. (|56|) . we find 

0(1 mm), (57) 



which is sufficient [23, 123 to evade current constraints on EP violation and fifth force. Note that 
the field redefinition Eq. H46(l can be taken as (pnew — ^/ci M In and the final result obtained 

above is still valid. In addition, let us stress here that the potential Eq. H49I) is different from those in 
Refs. [20I I21I I22*. '23] . Therefore it is interesting to further study this model more details and other aspects 
such as, the "thin-shell" effect. We will present these details in a separate paper 30]. 

The upshot of this subsection is that we illustrate briefly how the directly strong coupling between 
k-chameleon and matters (cold dark matters and baryons) is allowed while we cannot detect it from the 
tests of EP violations and fifth force searches on the earth or in the solar system today. The coupling 
constant (3 needs not to be tuned to an extremely small value and can be of order unity, in harmony with 
general expectations from string theory. The mass scale M can be of order Mpi. In other words, it needs 
not to be tuned like Eq. @ mentioned above as in Refs. IMillil. 



IV. COSMOLOGICAL EVOLUTION OF K-CHAMELEON AND THE COSMOLOGICAL 

COINCIDENCE PROBLEM 

In this section, we will study the cosmological evolution of k-chameleon. Just like the cases of 
quintessence |^ ^ and k-essence d in El 111, the key point is to find out its attractors during 
the evolution. During the radiation- and matter-dominated epochs, there are several sets of possible 
attractors in the k-chameleon model. Note that, as illustrated in Sec. IIII Cl the mass scale M can be of 
order Mpi. Thus, for simplicity, we take M — ^/SMpi from now on. 
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A. Radiation-dominated epoch 



In this epoch, the fraction energy density of matter flm — 0. One can see that aU terms describing 
the couphng between k-chameleon and matters can be neglected, and the master equations governing 
the system, namely Eqs. (I41ll - H44|l . reduce to corresponding ones for the ordinary k-essence case [TlLIT^. 
Thus, the cosmological evolution of k-chameleon in this epoch is completely the same as that of ordinary 
k-essence. A detailed study was presented in Ref. |12] for this case. Therefore we will not repeat here 
and only mention some key points and make some remarks. 




□ 



FIG. 1: A sample of function g{y), reproduced from |12|| . Here the new attractor C appears between yd < y < Vs 



Attractors 



During the radiation dominated epoch, three kinds of attractors, R (radiation), K (k-field), S (de 
Sitter) attractors (following notations of the k-essence model [Til H^) may exist. The location of the R 
attractor, i.e. yr = const., is determined by 



VrgiVr) = -3.9(yr), 



(58) 



which corresponds to yr < yd to ensure g > (positive pressure) , where y^ is the location of the function 
g{y) across the y coordinate axis, i.e. g{yd) — (see Fig. 1 of the present paper or Fig. 1 of Ref. ,12]'). 
The fraction energy density of k-chameleon is given by 



n 



(r) 



r^yr) = ~2g'{yr)yl 



(59) 



and the R attractor exists only if r'^{yr) < 1- The parameter of equation-of-state w^ijjr) = 1/3, and the 
k-chameleon mimics the radiation. Note that 57^'^'' = r'^{yr) has to be in the range 1 — 10% in order to 
satisfy the constraints from the big bang nucleosynthesis (BBN) [s^ . 
The location y^ of the second attractor K is determined by 

f^f^=r-2(y,) = l, (60) 

which implies that the k-chameleon dominates over other components. The parameter of equation-of-state 



-I 



2\/2 



1 



3 



const., 



and the scale factor is 



(61) 



(62) 
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If —g'^,y\l2 > 1, the solution describes a power-law inflation. Physically, this condition is equivalent to 
WtotiVk) = w^iuk) < -1/3. 

The third attractor S is defined by w^{ys) — —1 and fi^*' = f-^iUs) — 0. From Eq. 1^^ . we have 

g{ys)=9[ys, = 0. (63) 

Geometrically, S is a fixed point on the curve g{y)j which makes the tangent of g at ys passes through 
the origin (see Fig. 1). It is clear that ys always exists for the convex decreasing function g[y) unless 
ys oo. 

The combination of cosmologically relevant attractors during the radiation-dominated epoch can be 
one of three types: 

• R, S and no other attractors at < y < ys', 

• R, S, K plus possibly other attractors a,t y < yd] 

• R, S {no K attractor) and at least one additional attractor r(?) or k(?) (following notations of 
k-essence). 

The phase diagrams for these three cases have been drawn in Figs. 2, 3 and 4 of Ref. respectively. 
Obviously, one can see from these phase diagrams that the R attractor has the largest basin of attraction 
on the whole phase plane for all these three cases. 



2. Remarks 

The stability analysis of the attractors was done in Ref. j^. With a closer look, we find that the 
stability analysis is valid only for the R and D (dust) attractor (which may appear during the matter 
dominated epoch in the ordinary k-essence model). Actually, because the K and S are not determined 
by — Wm [following notations of Ref. [T^ . Wm denotes the state parameter of background matter 
(radiation or dust)], Eq. (24) of Ref. T^] is not vahd for them. In fact, Eq. (24) of Ref.]l3| should be 
changed to ddVl^/dN = 3w^{yk)Sil.^ and d5Q.^/dN ~ iSVl^ for K and S, respectively. Therefore, the 
stability condition for K attractor is ^/--g^jjfj2 > 2/3 while S is always unstable. This point can be seen 
clearly in the phase diagrams, i.e. Figs. 2-6 of Ref. jl^l: the phase flow can get near to S, but never 
reaches it, and the phase flow is then forced to leave it. 

The second remark is about the basin of attraction of the R attractor. In Ref. the statement "R 
attractor has the largest basin of attraction on the whole phase plane in the radiation-dominated epoch" 
was criticized by numerically analyzing two concrete models. In these two models the function D{X) in 
the Lagrangian has the following forms 

D{X) = -2.01 + 2VTTX + 3 X lO^^X^ - IQ-'^^X^, 

and 

D{X) = -2.05 -I- 2^J\ + f{X), 

where 

f{X) = X- lO^^X^ + IQ-^^X^ - IQ-^^X'^ + 10"2"X5 - 10"2'^XV2^ 

which arc first given in Refs. and respectively. (Note that in Refs. [llllT^[T3 |. the unit 3Mp; — 1 
was used.) However, we notice that the conclusion "R attractor has the largest basin of attraction on the 
whole phase plane in the radiation-dominated epoch" is based on another non-analyzable Lagrangian, 
i.e. Eq. (42) of Ref. [12, 

9{y) = 9giue{y) { 1 - ) , 

where g{y) is a function parametrized by five parameters: yr, g'r, yd, g'd and s^yd (see Sec. V of Ref. ^2 
for more detail). By choosing suitable parameters, as is shown in Sec. VA and VB of Ref. that the R 
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attractor has indeed the largest basin of attraction on the whole phase plane in the radiation-dominated 
epoch is possible. 

The third remark we would like to stress is about the implication of ys mentioned above. In the left 
side of ys, i.e. y < ys, one has > —1 and > while g" > 0. Because the g" is a continuous 
function, and the sound speed cannot diverge at ys, which implies that g" 7^ at y^- Therefore 
g" > should still hold in the right side of ys, i.e. y > ys- However, in the right side of y^, 5 ~ g'y < 0, 
one has < —1 and < 0. This is physically forbidden since the k-chameleon becomes unstable. As 
a result, the phase flow cannot pass across Any physically reasonable y must be less than ys, which 
can be seen clearly from Figs. 2-6 of Ref. [Tjl as well. Actually, a general discussion on this point has 
been made in Ref. which shows that a dynamical transition from the states with > — 1 to those 
with W0 < — 1 or vice versa is physically impossible. One can see that in fact, the condition y < ys is 
equivalent to > — 1 physically. 



B. Matter-dominated epoch 



In this epoch, the fraction energy density of radiation fir ~ 0, and Qm + ^(f> — 1 since the universe is 
spatially flat, as indicated by many astronomical observations |^. In this case, the terms due to the 
strong coupling between k-chameleon and matters cannot be neglected in the master equations Eqs. H41(l - 
H44|l . As will be seen shortly, due to the appearance of the coupling, the k-chameleon model will have 
a big difference from the ordinary k-essence model. In this subsection, we will find out all the possible 
attractor solutions and then study their stability. 



1. Attractors 



At first, we would like to point out that the D attractor [the scalar field tracks the dust {w^{yd) — 0) 
during the matter dominated epoch] and the S attractor, which may exist in the ordinary k-essence 
model [TllIT^ . are physically forbidden in our k-chameleon model, due to the existence of strong coupling 
between the chameleon and dust matters (cold dark matter and baryons). To see this point, let us note 
that in the cases of — 1, ^ct> — and Wtot — for the S attractor and Wtot = 'w<j> = 0, but 0^ ^ 
or 1 for the D attractor, the third term in Eq. H41|l is extremely large (for the case $7^ — > 0) while the 
first term in Eq. H42(l and Eq. (|43|l cannot be neglected. Thus, no solutions with y = yattractor = const. 
satisfying dy/dN — and d^tf,/dN = 0, dflm/dN = exist. Therefore the attractors D and S will no 
longer appear in the k-chameleon model. 

(k) 

However, we find that the k-chameleon dominated attractor K with fi^ — 1 can still occur in the 
k-chameleon model. Its characteristics is described by Eqs. (|60|I - H62|I . A remarkable feature we would like 
to stress here is that this solution describes a power-law inflation provided v^— 5^2/^/2 > 1. As we will see 
below, if this attractor is stable, this condition for accelerated expansion can be satisfied automatically. 

Except for the K attractor, we find that there is a new attractor solution in the k-chameleon model, 
which arises due to the strong coupling between the k-chameleon field and the dust matter. The new 
attractor is dubbed as C attractor (the letter C stands for "Chameleon" ) , and has some very interesting 
features. Let us describe the C attractor in some detail. 

Considering = -Kg' and K{(t>) = M"^ /cjp' [see Eqs. (jSHl) and and setting M = V^Mpi, from 
the Friedmann equation Eq. (|29|l one has 




Note that the relation a — f3/(t>, we have 



H=^^-. (64) 



(65) 



-g' 

which is independent of (/). Therefore it is possible to find an attractor solution C with y — yc — const. 
and fixed ri^''', il™'' = 1 — fi^''' and satisfying the master equations Eqs. H41|l - H44|l . 



2 
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and 



[^)jc + Sw4y,)n';;^^0. (67) 

Substituting Eqs. ^ and (|23 into Eq. we get 

V2/3 



^'f--i^V-9'. (68) 

which is determined only by the coupUng constant /3 and the function g{y) at y = Uc- Substituting fi^'^^ 
into Eq. H66|l . we have 

Comparing with r{yc) = Ms{yc) = VSMpi s{yc) and using Eq. (|38l) . we obtain 

2/?(2 + /3) 

9{yc)^9c^ r ■ (70) 

Geometrically, this means that the C attractor locates at the intersection of curve g{y) and the hyperbola 
h{y) = —2/3(2 + /3)/(9y) in the plot of g(ij) versus y. Because the asymptotes of the hyperbola h{y) 
are the two coordinate axes, and g{y) < when y > y^, this intersection always exists in the regime 
Uc > yd so that the k-chameleon contributes a negative pressure. On the other hand, note that the curve 
g{y) is monotonically decreasing while the hyperbola h{y) is monotonically increasing in the regime of 
Ud < y < ys, therefore there is only one intersection. In other words, the C attractor always exists and 
given a function g(jj), there is only one C attractor. 

Next let us have a look at the other physical features of the C attractor. From Eqs. (jfiSfl and H7()|l. the 
fraction energy density of the k-chameleon is 

r^(c) Qyc{-9c) .7-1 ^ 

> " 2(2 + /3)2- 

The usual restriction on fl^l^^ is < fl^j^^ < 1. However, as mentioned at the end of Sec. IIV A 21 the 
condition yc < ys has to be imposed. In that case, one has gc — g'cVc > while gc < 0, g'^ < and 
W(j,{yc) > — 1. By using Eqs. H68() and H70|l . we find that the bounds of Vt^^^ should be 

^ <^^:^<i- (72) 



2 + /3 



Note that (3 0(1), this result is quite interesting. For instance, the lower bound of rj^^-* are 1/3, 1/2, 
3/5 and 2/3 for /3 = 1, 2, 3 and 4, respectively. Therefore, in the C attractor, it is not strange that the 
fraction energy densities of k-chameleon and dust matters are comparable. In addition, from Eqs. H36|l . 
(EOJ), and {TTJ, one has 

MVc) = (73) 

Thus, the requirement w^{yc) > — 1 leads to the same lower bound to fi^'^'' as given in Eq. H72|) . Further- 
more, from Eqs. H44|l and H73|l . we find 

wtotiyc) = w^{yc)^^^^ = ^J^- (74) 



From the Einstein equation 

a 



Ptot(l + Swfot), (75) 
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the universe undergoes an accelerated expansion provided /9 > 1, which is equivalent to wtot{Vc) < —1/3. 
We can look at this from another angle. Because (j)c = \/2Xc = \/2/j/c is constant, one has (j) (xt. Then, 
from Eq. I|64|l . H oc t^^ and the scale factor a oc t^. If A > 1, the universe undergoes a power-law 
inflation. Let us find out the explicit expression of A for the C attractor. From Eq. H25() . one has 

PtotiVc) oc a-^^Ii+^-f^^)! = a-6/(^+2). (76) 

bmce n^^'^ and ni^;^ = 1 - nli;'^ are both fixed, p,f,{yc) and PmiUc) decrease in the same manner as 

PtotiVc) oc a"^/(''+2). Substituting a(xt^ and Eq. into Eq. we get A = (/?-|-2)/3 by comparing 
the power of t. In fact, from Eqs. (|64|) . H71|) one can find the same result 

a cx i('3+2)/3 (77) 

by using (p ~ (pet. In short, if /3 > 1, the C attractor solution describes a power-law inflation. 

In summary, in the k-chameleon model there may exist two attractor solutions, K and C, and no D 
and S attractor solutions in the matter-dominated epoch. If r'^{y) < 1 for any y < ys, the K attractor 
cannot exist. Thus, the combination of cosmologically relevant attractors during the matter- dominated 
epoch can be one of two types: 

• Only C and no K; 

• C and K. 



2. Stability analysis of the attractors 

As mentioned in the beginning of Sec. lIV A2l the stability analysis of the attractors K and C ought to 
be treated separately. We study the behavior of small deviations from the K and C attractor solutions 
one by one. 

(i) K attractor. In this case, fl^j^^ = 1, flm'' = 1 — fl^j^^ = 0. Substituting y{N) = yk + Sy and 
n^{N) = ^ + into Eqs. and hncarizing these equations, we get 



dSy^3w,iy,)~lf l 
dN 2 r'^ \^ ^ 2 ' 



iw^iyk) + L—^ 6n^. (78) 



Considering Eq. H3b|l. one has 'w^{yk) < since —g'f. > and < 0. Thus, the solutions of Sy and Sft^ 
decay only if 

3w4yk) + < 0. (79) 

Substituting Eq. H61() into it, the stability condition for the K attractor becomes 



^/-9',yl/^ > (80) 

Note that, if /? > 1 (the same requirement for the C attractor describes power-law inflation), Eq. (|80|l 
becomes y^—g'f,yl/2 > 1 which ensures the K attractor describes power-law inflation too. 

(ii) C attractor. In this case, ni^^ = 1- Qlj;^ Substituting y{N) = yc + 6y and fl^iN) = fl''^^ + SQ^ 
into Eqs. (|41|I - H43|I and linearizing these equations, we obtain 

^ ^ B,6y + B2Sn^, 

^^B,Sy + B,Sn^, (81) 
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where 



3 



i?4 = (2nl;^ - l) + \ ^ ^ V (82) 



From Eq. l|HT|l. we have 

^ - (5i + ^4)^ + (i?iS4 - B,B,)6n, = 0. (83) 
We can see that the solutions of Sy and Sil^ decay only if 

Bi+B4<0 and B1B4 - B2B3 > 0. (84) 

By using Eqs. (jSEJ, (EHl, (El, d)- IZSl and lIH3, and considering r{y) = Ms{y) = \/3Mpis{y), we can 
obatin 

4/3 + 2/?^ + 9(-gDyg 
(2 + /3)2 [4/3 + 2/32 + 9(-gi)y2] 



2 



i?3 = 



27g^'(-5i)y4 

3 [9(-ffD^yc + /3(/? + 2)g':] [2{I3 + 2)^ - 9(-g^)z/g 

4(2 + /?)4(-g^) 

/3[2(/3 + 2)^-9(-gOyg] 
6(2 + /3) (-5^)2/2 



B4 = 



From these quantities, we have 



[9(-g-)yg + 2/3(/3 + 2)] [2(/3 + 2)^ - 9(-ff-)yg] 
^154-^2^3- 4(2 + /3)V2/? • ^''^ 

Note that /? > 0, -g^. > 0, g'^ > [see Eq. (EHJ], and l^^''^ < 1, which implies 2{I3 + 2f -9{-g[)y'l > [see 
Eq. H71|l ]. Then it is easy to see that the stability conditions Eq. 184|l for the C attractor are satisfied. 
Thus, we conclude that the C attractor is stable. 
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C. Cosmological evolution of k-chameleon and the coincidence problem 



As mentioned in Sec. IIV Al in the radiation-dominated epoch, the behavior of k-chameleon field is 
the completely same as that of the ordinary k-essence without interaction between the scalar field and 
background matters a result, three kinds of attractors, namely, R, K and S, may exist. However, 

the radiation tracker, R, has been argued to has the largest basin of attraction on the whole phase plane 
so that most initial conditions join onto it and then makes this scenario became insensitive to initial 
conditions [U] ■ Yet the contribution of k-chameleon to the total energy density must not spoil the BBN 
or not dominate over the matter density at the end of the radiation-dominated epoch. It has been argued 
that if the contribution of the scalar field energy density in the R attractor satisfies ^3 

f^J'^ ~ IQ-^ - 10-\ (87) 

the scalar field (k-essence or k-chameleon) will not violate the constraint from the BBN. 

The new features appear in the matter-dominated epoch. As mentioned in Sec. IIV Bl in the matter- 
dominated epoch, the D and S attractors, which may occur in the ordinary k-essence model, are physically 
forbidden in the k-chameleon model, due to the strong coupling between k-chameleon and dust matters. 
Actually, in the k-chameleon model, two kinds of attractors may exist: one is the familiar K attractor 
and the other is a completely new attractor C. The new attractor C has some desirable features which 
can make the fraction energy densities of k-chameleon and dust matters (cold dark matter and baryon) 
be comparable. During the matter dominated epoch, the relevant attractors can appear in the following 
two possible sets: (1) only C and (2) C and K. If /? > 1, the universe will undergo a power-law inflation, 
regardless the k-chameleon enters into the C or K attractor. Suppose that during the radiation dominated 
epoch, the universe enters into the R attractor, after the onset of matter domination, it will enter into 

(r) 

the C attractor with a large possibility since the values of fl^ and are required to be somewhat 
small in order to satisfy the BBN constraint. In particular, one can adjust the model so that the K 
attractor does not exist (for instance if r'^{y) < 1 ior y < i/s), thus the universe has to enter into the 
unique C attractor where the fraction energy densities of the k-chameleon and the matters Qm are 
fixed and keep comparable forever. In this way the k-chameleon model leads to a natural solution to the 
cosmological coincidence problem. 

The evolution of the k-chameleon heavily depends on the function g{y) and other components of the 
universe like radiation and dust matter. In Fig. 2 we plot a sketch of possible phase diagrams of the 
evolution of the k-chameleon, where only two attractors R and C appear during the evolution of the 
universe. In this plot, we expect that during the radiation domination epoch, for most initial conditions, 
the k-chameleon is attracted to the R attractor satisfying the constraint H87(l . In this epoch, the k- 
chameleon mimics the equation of state of the radiation component of the universe. With the increase of 
flm for the dust matter component, the k-chameleon will no longer track the radiation component due 
to the interaction between the k-chameleon and dust matter. During the matter domination epoch, the 
fl^ will continue to decrease until y reaches to yc, which can be seen from (|42(l and l|43(l . Beyond ?/c, ^<p 
will increase and y decreases toward to yc and passes through it, which can be seen from 1)41(1 . When 
y decreases to some value (< yd), it increases towards to y^ again. After several such processes, finally 
the k-chameleon is expected to reach the stable attractor C, where the fraction energy densities of dust 
matter and dark energies are comparable and the universe undergoes an accelerated expansion. 

In order to have a better picture of the k-chameleon model, it is helpful to do some numerical analysis. 
The astronomical observations, such as SNe la 0,|3,|3|, WMAP 0|, suggest that fi^ ~ 0.7, fi^ — 0.3, 
and < —0.76 at 95% C.L. today. If adopt ^l^{yc) = 0.7 and impose the constraint w^{yc) < —0.75, 
we see from Eqs. JZU and that 2.2 < /3 < 4.7 has to be obeyed. From Eqs. and we find 
that w^{yc) ~ -0.86 and wtotiyc) = -3/5 for /3 = 3, while Wcj,{yc) ^ -0.95 and Wtot{yc) = -2/3 for 
(3 = 4:. In these cases, the universe undergoes an accelerated expansion as a cx t^^"^ and for /3 = 3 and 
4, respectively. 

In the C attractor, the kinetic energy term Xc = l/^c fixed at a constant value. It is possible to 
design the function g{y) to get a somewhat large yc- For instance, if Q^{yc) = 0.7 and /3 = 3, we have 
from Eq. lO that y^ ~ 100 and ~ 10"^ while ~ -4 x 10"^ (the unit SM^, = 1 has been used 
here). Comparing g'^ ~ —4 x 10^"* with a particular example g'^ ~ —5 x 10^"^ in Ref. can see 

that this value is reasonable, since yc — 100 ^ yd = 17. Therefore, in this example, the k-chameleon 
indeed can be treated as a canonical chameleon approximately. As a result, we cannot detect it from the 
tests of EP violations and fifth force searches on the earth or in the solar system today, although it is 
strongly coupled to ambient matters, as illustrated in Sec. IIII Cl 
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FIG. 2: A sketch of possible phase diagrams for the case where only the R and C attractors appear. 

Note that in order to exclude the K attractor, one has to adjust the model so that r^(y) < 1 and 
decreases monotonically in the region yr < y < Ds p^ . On the other hand, one requires that the C 
attractor exists in the region yd < y < ys- One may wonder whether or not these two conditions can 
be met simultaneously. To see this, let us take an example. We have from Eq. H66|l that r^(yc) = 
[(2 + I3)n'^^^ - PY/i'in^f). If /3 = 4 and n'^^^ = 0.7, one then has r^{yc) ^ 1.4%. In the radiation 
dominated epoch, in order to satisfy the constraint from the BBN, r'^{yc) — ^If^'' should be in the region 
1 — 10%. So we see that those two conditions can be satisfied, if r'^{yr) > 1-4% and the decreasing of 
r^(y) is sufficiently slow so that t/c and ys can be somewhat large values. Note that yc is always less than 
ys since ys locates at r'^{ys) = 0. 



V. CONCLUSION 

Recently a chameleon mechanism has been suggested 0, 0, in which a scalar field 

(chameleon) can be strongly coupled to ambient matters, but it still satisfies the constraints from the 
fifth force and EP violation experiments on the earth and in the solar system. In this paper we have 
combined the chameleon mechanism to the k-essence model of dark energy and have presented a k- 
chameleon model. During the radiation dominated epoch, the evolution of the k-chameleon is the same 
as that of the ordinary k-essence, and three kinds of attractors, R, K and S, may appear. One can 
construct a model where the R attractor has the biggest basin of attraction so that for most initial 
conditions, the universe will be attracted to the R attractor. During the matter dominated epoch, the D 
and S attractors, which may appear in the ordinary k-essence model, are forbidden in the k-chameleon 
model, due to the strong coupling between the k-chameleon and background matters (cold dark matter 
and baryons). Except for the familiar K attractor, a new attractor, dubbed C attractor, exists in the 
k-chameleon model. In the C attractor, the fraction energy densities of the chameleon field (dark energy) 
and the dust matter (cold dark matters and baryons) are fixed and comparable, and the universe enters 
into an accelerated expansion phase in the power-law manner if the coupling constant (3 > 1. We can 
adjust the model so that the K attractor does not exist, the universe then has to enter into the C attractor. 
Thus the k-chameleon model provides a natural solution to the cosmological coincidence problem. 
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